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ABSTRACT 

Context. It is widely accepted that solar flares are manifestations of magnetic reconnection events taking place in the 
solar atmosphere. Several aspects of these events remain unclear, although many efi'orts have been devoted to the 
investigation of magnetic field configurations at fiare occurrence sites. 

Aims. In this work, we have studied the temporal evolution of some properties of a sample of superactive regions with 
the aim to single out the most significant for flare activity forecasting. 

Methods. We have investigated properties of 14 superactive regions, observed between January 1st 2000 and December 
31st 2006 with MDI/SOHO instrument and characterized by a particularly intense fiare activity during their passage on 
the solar disk. We have analyzed the temporal evolution of fractal and multifractal properties of photospheric magnetic 
fields, namely the generalized fractal dimension and the cntribution and dimensionality diversities, which describe 
geometrical properties of the magnetic field, as well as the potential unstable volumes of magnetic discontinuities 
above the studied ARs, which may provide information about the magnetic field configuration in upper layers of the 
atmosphere. Correlations of these quantities with the fiare index, which provides information about the flare activity 
of a region, have also been estimated. 

Results. We found that in 50 % of our sample the generalized fractal dimension is correlated with the fiare index 
computed over windows of 50 hours, while the contribution diversity and the dimensional diversity are anticorrelated 
with the same index. A clear increase of the potential unstable volume of magnetic discontinuities in the corona is 
observed before the phases characterized by more frequent and intense fiares. We also found that the free energy 
distribution functions of unstable volumes of the analyzed superactive regions can be fitted with straight lines whose 
slope is larger than the values found in previous works for less active magnetic regions. 

Conclusions. The generalized fractal dimension and the potential unstable volume of magnetic discontinuities are the 
most suitable for statistical investigations of relations with flare activity over longer (50 hours) and shorter (few hours) 
time intervals, respectively. 

Key words. Sun: activity - Sun: magnetic flelds - Sun: fiares 



^ 1. Introduction Many statistical studies have been devoted to the in- 
vestigation of physical properties and topology of the 

It is largely accepted that solar flares are magnetic recon- magnetic field in the photosphe re, compatible with the 

J5 . nection phenomena that can occur at different scales of time events (see iSchriTver \\2mA [2&59L and references t h erein) , 

■ - • and energy (|Aschwanden et al. 1 120001). Although the flare but the results are contro versial (|Georgoulis I [2001 . 

intensity and frequency are proportional to the magnetic i Rpmano fc Zuccarello I (fMl (Paper I, hereafter) investi- 

ficld strength involved in these events, it is still unclear g^ted the correlation of flare productivity and some prop- 

what the magnetic configurations are that drive these en- gj-ties of 26 superactive regions. They found that the devia- 

ergy releases. In the past decade the theory and the nu- tions from magnetic polarity balance and the length of the 

merical simulations of magnetic reconnection have shown inversion lines characterized by a horizontal magnetic gra- 

that the dissipation of magnetic energy, responsible for the dient greater than 0.3 G km"! seem to have a major role in 

flare activity occurs along unstable magnetic discontinu- producing a magnetic field topolo gy able to trigger flares, 

ities ([Parker I [19881). The formation and disappearance of This result has been confirmed bv FSong et al"[ ^M), who 



these discontinuities, also named current sheets, is a con- analyzed 174 active regions and employed a ordinal logis 

tinuous process which occurs in all magnetic concentrations tic regression model to predict sola r flares. Different con 

as in the solar active regions. Due to the thickness of these elusions were drawn by iLeka fc Barnes I (|2007D . who ap- 

discontinuities, below current telescope resolution, many p^jed discriminant analyses to almost thirty properties of 

authors have tried to deduce their presence and location the photospheric magnetic field estimated for 496 active 

from the global topology of active regions. regions. Their results indicate that available data do not 

provide a basis for reliable flare prediction. However, these 
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by including descriptors of the self-similar properties of the 
photospheric magnetic field in the analyses. 

Self-similar properties of the solar magnetic field are 
a signature of the turbulent physical processes which 
determine the generation, evolution and dissipa tion of 
magnetic structures (see for instance ICattanc o et al 



20031 lAbramenko fc Longcopel 120051 : IStein fc Nordlund 



20061 and references therein) and have b een investi- 
gated in several studies (see for instance iGeorgouIisI 
l2005l for a review). In this context, investiga- 
tion of fractal and multifractal properties of photo- 
spheric magnetic features provide important informa- 
tion about the ph y sical prope r ties o f convective plasma 
(iBalke et all 119931: iTao et al. I 1X9951: iNesme-Ribes et al l 



19961). of the magnetic field JStenflo fc Holzreuter 1 120021: 



Pietarila Graham et al. |[2009l ) and about the physical pro- 



cesse s that govern their interaction (iLawrence et al. Ill995l 
119991: iBushbv fc Houghton I l2005t ICrouch et arTl2007D . 

Studies carried out on short temporal ranges (hours-days), 
have showed that fractal and multifractal properties vary 
with the evolution of the photospheric magnetic field of 
the regions analyzed (e.g. iGeorgouhs 1 120051; [Abramenkol 
I2005at IConlon et alTI 120081: iHewett et al. 1 120081) . Several 
authors have also investigated the self-similar properties 
of AR s in connectio n with their flare productivity. For in- 
stance [Abramento] ([2005b) found that the magnetic power 
spectrum of 16 investigated ARs b ecomes steeper with 
the increase of the flare productivity. lAbramenko I (|2005al ) 
also compared the multifractal properties of two regions 
characterized by different flare activity levels and found 
that the more flare productive region showed a higher 
level of multifractah t y. A similar investigation was carried 
out by iGeorgouiisI ( 20051 ). who also presented the mul- 
tifractal analyses as more suitable than the fractal one 
for flare forecasting. Statistical analysis of large samples 
of data have shown th at variations of fractal properties 
of ARs in photosphere (iMeunier I l2004f ) and low chromo- 
sphere (jCriscuoli et al. I [2007f) do not present clear trends 
with t he activity cycle. On the contrary, McAtcer et al. I 
(|2005|) reported a correlation between variations of photo- 
spheric fractal properties of ARs with flare productivity. 

Precursors to flares have also been looked for in the so- 
lar corona, where the sites of magnetic reconnection are 
usually located. Variations of magnetic field characteris- 
tics and topology in connection with oc currence of eruptive 
phen o mena have been investigated (e.g. lUgarte-Urra et al. I 
2007t ISchrii ver 1 12009|) . Among the most recent studies, 
Barnes fc Leka I (j2006l) applied discriminant analyses to 



several properties of coronal fields above three active re- 
gions and concluded that coronal field descriptors are more 
efficient than photospheric ones to distinguish between 
flaring and non-flaring regions. Self-similarit y is also a 
chara c teristic of coroiial mag netic fields (e.g. IGeorgouiisI 



20051: [Belanger et al. 2007r and references therein). 



Aschwanden fc Aschwandenl ()2008al lbl) measured the frac 



tal dimension of flare arcades and foun d it to b e linearly 
correl ated with the flare flux in EUV. lAbrame nko et al. 
()2008l ) estimated the intermittency of the magnetic field 
of an AR and found a time delay of about one day be- 
tween the evolution of intermittency i n photosphere and 
in corona. IVlahos fc Georgoulis I ()2004[ ) investigated prop- 
erties of the magnetic discontinuities of the solar corona, 
i.e. regions characterized by magnetic field variations oc- 
curring on short scale lengths and where the nonpotential 



magnetic energy is accumulated. While most of these dis- 
continuities remain relatively stable, a fraction of them may 
become unstable, thus causing rapid magnetic reconnection 
that releases the free magneti c energy and provides the en- 
ergy content for flaring events. IVlahos fc Georgoulisi (|2004l ) 
studied an active region characterized by no flare or subflare 
occurrence and an active region showing a number of sub- 
flares but no major flares over the course of several hours. 
They used linear force-free field extrapolations to calcu- 
late the 3D magnetic configuration and various selection 
criteria to identify these discontinuities. They found that 
the free magnetic energies and volumes contained in these 
discontinuities obey power-law distribution functions, with 
power-law indices for the free energies in the range [-1.6, 
-1.35]. 

In this paper, we extend the work presented in Paper I 
by investigating other properties of photospheric and coro- 
nal magnetic fields of the same sample of superactive re- 
gions. We consider the temporal evolution of some frac- 
tal and multifractal properties of the photospheric field, 
namely the generalized, the cont ribution and the diversity 
dimensions (jConlon et al. I [20081 ) (see §3.1 for the defini- 
tion of these parameters), estimated by the analyses of MDI 
magnetograms. We also investigate the evolution of the free 
energy distribution and the potentially unstable volume of 
the magnetic discontinuities (Vlahos & Georgoulis 200^ 
of the coronal field (see §3.2) as derived by force-free field 
extrapolation. We investigate the temporal correlation with 
the flare activity in order to single out the indicators that 
best describe the flare productivity potential of a region. 
We choose these indicators because they describe the topo- 
logical complexity of the magnetic field in different layers 
of the solar atmosphere. 

This paper is organized as follows. The data-set ana- 
lyzed and the indicators employed for this study are de- 
scribed in §2 and §3, respectively. In §4 we present the re- 
sults we obtained and in §5 we discuss them and draw our 
conclusions. 



2. Data 

According to criteria explained in Paper I, we selected 26 
superactive regions (ARs) which appeared on the solar disk 
during the middle of the 23th solar cycle, between January 
1st 2000 a nd December 3 1st 2006 and characterized by a 
flare index (|Li et al. Il2004[ ) greater than 500, computed over 
time intervals of 50 hours. The X-ray flare index I computed 
over the time interval At is defined as: 

/(At) = O.l Bh + ^a + 10 Mj + 

h i j 

+ 100 X ^ Xk (1) 

k 

where Bh, Ci, Mj and Xk refer to the soft X-ray classifica- 
tion based on the flu x in the 1 -8 A range of the spectrum 
(Tandberg-Hanssen fc Emsliel il988). Depending on the X- 
ray flux, flares are classified as B, C, M or X, where B are 
the less energetic flares (lO'^erg • cm~^ and X are the 
most energetic (lO^^erg • cm^^ ■ s^^). The sums in eqnllj 
run over the number of flares of that class produced in an 
active region in the time interval At. 

For each region we used MDI/SOHO line-of-sight full- 
disk magnetograms taken during the time interval between 
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Table 1. Superactive regions selected, correlation coefRcients between the investigated indicators and the flare index, 
average value of indicator as- 



NOAA 


First Flare 


Last Flare 


N mag 


Ds 


Da 


Cd 




< as > 


9026* 


Jun. 1,2000 01:01 UT 


Jun. 10,2000 178:19UT 


48 


0.68 


-0.62 


-0.34 


-0.38 


-1.6 ± 


0.2 


9236 


Nov. 22,2000 04:59 UT 


Nov. 27,2000 00:21UT 


82 


0.75 


-0.68 


-0.59 


0.54 


-1.7 ± 


0.3 


9690 


Nov. 6, 2001 08:34 UT 


Nov. 15, 2001 17:29UT 


84 


0.66 


-0.77 


-0.65 


0.55 


-1.7 ± 


0.4 


9733 


Dec, 8, 2001 02:56 UT 


Dec. 19, 2001 01:24 UT 


87 


-0.72 


0.62 


0.1 


0.44 


-1.55± 


0.4 


10030 


Jul. 10, 2002 01:56 UT 


Jul. 20, 2002 12:44 UT 


94 


-0.02 


0.04 


0.29 


-0.40 


-1.8 ± 


0.4 


10039* 


Jul. 23, 2002 00:47 UT 


Aug. 3, 2002 19:11 UT 


33 


-0.02 


0.54 


0.66 


0.04 


-1.7 ± 


0.2 


10314 


Mar. 15, 2003 03:44 UT 


Mar. 22, 2003 23:03 UT 


46 


0.23 


-0.57 


-0.38 


0.04 


-1.7 ± 


0.2 


10365* 


May 25, 2003 20:38 UT 


Jun. 2, 2003 16:01 UT 


68 


0.81 


-0.78 


-0.66 


-0.29 


-1.7 ± 


0.3 


10375* 


Jun. 1, 2003 04:00 UT 


Jun. 15, 2003 12:35 UT 


83 


0.5 


-0.03 


0.14 


0.46 


-1.65± 


0.3 


10486* 


Oct. 22, 2003 05:21 UT 


Nov. 5, 2003 10:56 UT 


80 


-0.2 


-0.03 


0.42 


-0.93 


-1.45± 


0.2 


10656* 


Aug. 6, 2004 23:01 UT 


Aug. 19, 2004 14:16 UT 


76 


0.88 


-0.66 


-0.59 


0.84 


-1.7 ± 


0.3 


10696 


Nov. 2, 2004 16:42 UT 


Nov. 11, 2004 16:02 UT 


91 


-0.11 


0.68 


0.54 


0.22 


-1.7 ± 


0.2 


10720 


Jan. 11, 2005 20:46 UT 


Jan. 23, 2005 05:23 UT 


81 


-0.03 


0.1 


0.06 


0.38 


-1.8 ± 


0.3 


10786* 


Jul. 1, 2005 05:07 UT 


Jul. 15, 2005 17:15 UT 


74 


0.5 


-0.35 


-0.41 


-0.77 


-1.8 ± 


0.2 



the first and the last flare, with a time resolution of 96 
minutes and a plate scale of 1.98 arcsec/pixel. We con- 
sidered only those magnetograms where the longitudinal 
angular distance of the active region from disk center was 
less than 41 degrees. Among these 26 ARs we neglected 
those fo r which MDI sho wed a saturation effect in sunspot 
umbra (|Liu et al. Il2007f) . thus reducing the analysis to 14 
ARs. For each region we considered a sub-field centered on 
the baricenter of the AR such that the maximum magnetic 
flux imbalance was less than 20%. Sub fields of size smaller 
than 100 X 100 pixels were discarded. The number of magne- 
tograms for each AR selected according to these criteria is 
given in Table [TJ Geometrical projection was comp ensated 
for b y the ideproject Solar Soft routine ( Freeland et al. I 
fl998l) . 

For 7 ARs (marked with a in Tab. [T]), which have 
been observed by TRACE, we also considered high resolu- 
tion images at 171 A or 195 A in order to investigate their 
coronal configuration. 



3. Data analysis 

From each line-of-sight magnetogram sub-field we com- 
puted 5 parameters, three of which describe fractal and 
multifractal properties of the photospheric magnetic field 
and two which describe the behavior of the upper level of 
the solar atmosphere above the selected ARs. Then we in- 
vestigated the correlation between these properties and the 
flare index computed over temporal windows of 50 hours. 



3.1. Fractal and multifractal indicators 

To investigate the fractal properties of ARs, we estimated 
the generalized fractal dimension by using the box- 
counting technique. The computation of the Contribution 
Diversity Cd and Dimensional Diversity Dd allowed us to 
investigate multifractal properties. As already pointed out 
in §1, the literature concerning fractal and multifractal 
studies in solar physics is quite extensive; here we sum- 
marize the main points. 

The most common method of estimating the fractal di- 
mension of an object given its image is box counting. It 



consists of covering the image with boxes of different sizes 
e and then estimating the slope of the linear relation 



IniV(e) ^ DVanln - + C 



(2) 



where N{e) is the number of boxes that cover the object, 
C is a constant and D is the fractal dimension. More in 
general, given a measure P, the generalized dimension is 
defined as 



1 ln/g(£) 
Dq = lim —-^ 

q-l c^o Ine 

where g is a real number, 



(3) 



(4) 



and we have introduced the more compact notation iV^ in- 
stead of N{e). As Pi{e) we consider the normalized mag- 

netic fiux — — , where j runs in a box of size e and (ptot 
is the total flux in the image. The variation of the exponent 
q allows us to investigate different spatial scales: negative 
values of q enhance the smallest spatial scales whereas posi- 
tive values enhance the largest spatial scales. Dq is a mono- 
tonic decreasing function of q. The Renyi dimension Dq is 
the fractal dimension. Instead of Z?o, we investigated the 
temporal variation of the generalized dimension Dg, where, 
according to equation [31 



Da = — lim 



In /8(e) 



7 e^o Ine 



(5) 



We choose D^ since it is less sensitive to spatial resolu- 
tion and pixelization effects (|Lawrence et al. I[T99^ . 
The strength of the measure is defined as 



Ine 



Q!o = lim 



(6) 



where r;,,,(e) - Y^^l^^l^^y ■ 

The multifractal spectrum F{aq) is defined as 



Fiaq) = lim 

^ e-+o Ine 



(7) 
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It can be shown that, under the assumption F"{aq) < 
<^ dDg/dq < 0, 



F{aq) = qoq - {q- l)Dg 



(8) 



aq is the exponent which describes the scahng of the 
measure at a cer tain position of the domain, while F{aq) is 
its distribution (jHalsev et al. I [19861 ). For a bi-dimensional 
image < < 2, where = in the case of a sharp spike 
and aq — 2 for a homogeneous image. Therefore, low values 
of Uq indicate denser regions, while high values indicate 
more rarefied regions, and the value of aq describes th e 
degree of 'spikiness' of an image (|Chappell fc Scalol[200ll) . 
The set of points of the domain with the same value of 
aq has fractal dimension F{aq). The spectrum is a concave 
function of aq and its maximum is the fractal dimension D^. 
For a single fractal of fractal dimension aq — F{aq) = 
Df,yq. The c ontribution diversity Cd and the dimensional 
diversity Dd (jConlon et al. 1120081 ) are defined respectively 



Therefore, the contribution diversity and the dimen- 
sional diversity describe the variety of the concentration 
of the different structures in an image (spikes, filaments, 
clumpy structures) a nd of the their complexity. Following 
IConlon et al. I (|2008D . the analysis was restricted to g > 0, 
since numerical errors are larger for negative values of the 
exponent q. 

3.2. Volume and free magnetic energy of unstable magnetic 
discontinuities 

For each magnetogram we performed a linear force-free 
field extrapolation to determine the 3-D magnetic field 
configuration. We applied an extrapolation cod e based on 
the method introduced bv I Alissandrakis 1 (jl98l[ ). which as- 
sumes that the magnetic field is force-free not only in the 
corona but also in the lower levels and that it vanishes at 
infinity. In order to minimize the oscillatory behavior of the 
Fourier solutions, we put the original arrays of the data for 
the boundary in bigger arrays of double dimensions, be- 
fore applying the extrapolation procedure. A consideration 
of the size of the integration box gives a resonant value of 
the force-free parameter. Below this value the periodicity 
effects inherent to the Alissandrakis method are reduced. 
Therefore, for each AR we fixed the force-free parameter 
just below the resonant value for the magnetogram with 
the largest field of view among the magnetograms of the 
same AR. Moreover, in order to investigate the sensitivity 
of our results to variation of this parameter, for 7 AR we 
also performed extrapolatons using the force-free parame- 
ters which best fit the coronal loops observed at 171 A or 
195 Awith TRACE. 

Fo llowing the method used by IVlah os fc Georgoulis I 
(120041 ) we applied a selection criterion based on the Parker 
angle in order to identify potentially unstable regions in- 
side each AR. In particular, we considered as potentially 
unstable the location where the angular difference between 
two adjacent magnetic field vectors was greater than 30°. 
Then, we computed the normalized total volume occupied 
by all the potentially unstable regions: 



where N is the number of the potentially unstable regions, 
rii is the number of adjacent locations inside each region, 
A is the pixel size of the magnetogram and Sh is the height 
step of the force-free field extrapolation. The normalization 
to the total volume of the box It in which the extrapolation 
was performed takes into account the variable sizes of the 
sub-field magnetograms analyzed. 

We also measured the free magnetic energy Ei in each 
unstable region i consisting of a set of neighbouring grid 
points j where the instability criterion is met: 



E, = 



X^Sh 
"8^ 



2 



B 



(11) 



with Bffj and -Bp/j the magnetic field values deduced 
by the linear force-free field and the potential field, re- 
spectively. Note that this f ormula, which has been widely 
used in the literature (e.g . iJing et aT~l l2009t iBobra et al. I 
120081: iKusano et al. 112002), diff ers from the one reported 



Vlahos fc Georgoulis 



( 2004 ) for reasons not known to 
these authors. 

Considering the distribution function of the total free 
energy Ei in the selected discontinuities, we obtained a 
power law with a different power law index aE for each 
magnetogram. This indicator was estimated by a linear fit 
to the logarithm of the total free energy distribution. Since 
the shapes of the distribution functions evolve with time 
(see for instance Figs. [H [H and we employed an au- 
tomatic technique (jMcAteer et al. |[2005l) to determine the 
linearity ranges and perform the linear fits. 

The difference between 1^ and a e estimated using force- 
free parameters obtained by the best fit of the observed 
loops and using the resonant values are on average 10% 
for both quantities. These differences are lower than tem- 
poral variations measured for /„ and are within the error 
associated with the estimation of as (see §4.1), thus show- 
ing that the choice of using as a force-free parameter the 
magnetogram resonant value produces reasonable results. 



4. Results 

4.1. Temporal evolution of the indicators 

We studied the temporal evolution of Dg, Cd, Dd, ly and 
aE for each superactive region. For some indicators, we 
found these variations to be dependent on the evolution of 
the ARs. 

We found that for all the active regions characterized 
by an emerging magnetic flux, the multifractal spectrum 
'contracts' and gradually shifts toward higher values of 
both aq and F{aq). As shown for instance in Fig. [T] for 
NO A A 9236, the MDI magnetograms taken on 22, 23 and 
26 November 2000 (first row) show the emergence of flux 
in photosphere, and the corresponding plots of multifractal 
spectra (second row) show that the values of both F{aq) 
and aq increase, while the ranges of values spanned by 
th ese quantities gradua lly reduce. As already pointed out 



bv lConlon et al. I (|2008f ). who found similar results, this re- 



(10) 



fleets the fact that magnetic flux typically emerges frag- 
mented in small-size features, which coalesce and evolve 
into larger and more complex structures. On the contrary, 
during the decay phase both aq and F{aq) decrease while 
the spectrum broadens, as illustrated in Fig. [2] for NO A A 
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Fig. 1. Top: Details of MDI magnetograms, acquired on 22 November 2000 04:51 UT (left), 23 November 2000 11:20 
UT (center), 26 November 2000 19:11 UT (right), showing the emergence phase of NOAA 9236. Middle: corresponding 
multifractal spectra. Bottom: corresponding distribution functions of total free energy. 




0,5 1.0 !,5 2-C 0,5 1,0 1,5 2,0 0,5 1,0 1,5 2,0 




25,0 25.5 26,0 26,5 27,0 27,5 28,0 25,0 25,5 26.0 26.5 27.0 27.5 28.0 25.0 25.5 26.0 26.5 27.0 27.5 28.0 

Total free energy (erg) Total free energy (erg) Total free energy (erg) 

Fig. 2. Top: Details of MDI magnetograms, acquired on 05 June 2000 00:00 UT (left), 08 June 2000 04:48 UT (center), 
10 June 2000 16:00 UT (right), showing the decay phase of NOAA 9026. Middle: corresponding multifractal spectra. 
Bottom: corresponding distribution functions of total free energy. 



9026. This probably reflects the fragmentation of the mag- 
netic field into less and less complex small-scale regions as 
confirmed by the magnetograms shown in the first raw of 
Fig.l 

For ARs in the emerging phase the fractal dimension 
Ds increases with time, while both multifractal indices Dd 
and Cd decrease, as shown for instance in Fig[4] (a). The 
opposite is found for NOAA 9026 during its decay phase. 



as shown in FiglH (b). In these plots t=0 corresponds to the 
time when the first flare was registered. 

For regions observed around the maximum of their life, 
the evolution of the multifractal spectrum is more com- 
plex. In general, we found that the fractal dimension stays 
fairly constant, while larger variations of multifractal indi- 
cators are observed at the emergence of new magnetic flux. 
This is particularly evident, for instance, for NOAA 10030, 
as shown in the appendix. Some cases also show that at 
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25.0 25-5 26.0 26-5 27.0 27,5 28, C 
Total free energy (erg) 





25-0 25,5 26-0 26,5 27.0 27,5 28,0 
Total free energy (erg) 





25-0 25,5 26-0 26,5 27,0 27,5 28,0 
Total free energy (erg) 



Fig. 3. Top: Details of MDI magnetograms, acquired on 11 December 2001 06:27 UT (left), 14 December 2001 09:36 UT 
(center), 16 December 2001 09:36 UT (right), showing the temporal evolution of NOAA 9733, Middle: corresponding 
multifractal spectra. Bottom: corresponding distribution functions of total free energy. 




0,0020 
' 0,0015 
0,0010 
0,0005 









(a) 





150 200 
Time (hours) 



(b) 



Fig. 4. Temporal variations of the estimated parameters for (a) region NOAA 9236 and (b) NOAA 9026 (b). The squares 
in the last panel of (b) indicate the /„ values deduced from the force-free field extrapolations performed using a best fit 
criterion with the loops observed by TRACE at that time. 



the very beginning of the decay phase, that is, during the 
hours after the maximum of flux is measured, the fractal di- 
mension continues to increase and multifractal indicators to 
decrease. Figure |3] displays for instance the magnetograms 
and the corresponding multifractal spectra of NOAA 9733 
at the beginning of the observations (left), at the moment 
in which the maximum total flux was measured (center), 
and at the end of the observations (right). Comparison of 
images in the left and middle columns shows that with the 
emergence of new magnetic flux, both the number of com- 



pact larger features and of fragmented smaller features in- 
creases. The former result from the coalescence of 'old' flux 
into larger and compact features, whereas the latter are 
the signature of 'new' emerged flux. As a consequence, the 
multifractal spectrum broadens and /(amm), i-e. the frac- 
tal dimension of the largest scale structures, decreases. The 
right column shows that, at the very beginning of the de- 
cay phase, the small features coalesce into larger and more 
complex regions, so that the spectrum 'contracts' and shifts 
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toward higher values of both f{aq) and ag, although the 
total amount of magnetic flux is decreasing. 

On the other hand, the temporal evolution of /„ does 
not seem to be strictly correlated with the evolution of 
the magnetic flux in the photosphere. Although in some 
cases during the emerging phases increases, as shown in 
FigHl (a), and during the decaying phases it decreases, as 
shown in Figdl^b), in other cases the amount of magnetic 
discontinuities changes independently from the magnetic 
flux observed in the photosphere. This result confirms that 
the coronal field extrapolated using the linear force- free ap- 
proximation provides some independent information from 
the magnetic field above the level of the measurements. 
We also found that for most of the regions analyzed, ly 
does not vary monotonically, but it is rather a discontinu- 
ous function of time. Square symbols in the plot of in 
FigHl (b) represent the total unstable volume computed by 
using a force-free parameter derived from the comparison 
with TRACE data. The values obtained in this way are 
very close to the ones obtained when the force-free field 
parameter is fixed just below its resonant value. 

The analysis of the temporal evolution of the distribu- 
tion functions of total free energy of unstable volumes shows 
that during the emergence phase the peak of the distribu- 
tion, as well as the linearity range, shift gradually toward 
higher energy values, as shown for instance in FiglT] and 
FigO On the contrary, temporal variations of the energy 
distribution of region NOAA 9026 (Figl2]), observed during 
the decay phase, do not seem to have a particular trend 
with the evolution of the region. Variations with time of 
the slope of the linear fits to the distributions [ue indica- 
tor) are within the error bars and are not correlated with 
the evolution of the regions. 

4.2. Correlation with the flare index 

In order to investigate the correlation of the magnetic field 
topology with the flare occurrence we reported the flare 
index and the other parameters in the same plots for each 
superactive region (see for instance Figs. [J] (a) and (b) and 
Appendix A). We represented each single flare with a black 
bar whose height is proportional to the X-ray flare peak 
observed by GOES and the total of the flare index over 50 
hours with grey bars. We used the logarithmic scale to show 
the magnitude of different flares in the same plot. 

We estimated for each AR the correlation coefficients 
between the considered indicators and the flare index com- 
puted over windows of 50 h. Results are reported in Table 
[T] (columns from 5 to 8). 

In approximately 50% of the ARs we found that Dg 
is positively correlated with the flare index: 7 ARs have 
a correlation coefficient larger than 0.5, with the highest 
value of 0.88 for NOAA 10656. For the same ARs (except 
for NOAA 10375) we obtained that Dd and Cd are an- 
ticorrelated with the flare index. The smallest correlation 
coefficients between the multifractal parameters and the 
flare index are -0.78 and -0.66 for NOAA 10365. For other 
regions the correlation between fractal and multifractal in- 
dicators is almost null. A comparison to total flux variation 
(see Appendix A) reveals that most of these regions have 
been observed at about the maximum of their evolution. 
The peculiar values found for NOAA 9733 are due to the 
large decrease of the flare index at the beginning of the de- 
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Fig. 5. Maximum of the flare index versus average gener- 
alized fractal dimension Dg for the 14 ARs analyzed. 



cay phase, during which, as already pointed out in previous 
paragraph, increases while Dd and Cd decrease. 

Figure[5]displayes the relation between the average frac- 
tal dimension Dg and the maximum measu red for the flare 
index. In agreement with previous works (jMcAteer et al. I 
120051: 1 Abramenko Il2005bl ). the plot suggests that more com- 
plex regions have a higher probability of developing more 
intense and a larger number of flares. Nevertheless, the large 
scatter of data does not allow us to derive an analytical re- 
lation between Dg, and the flare index. Note that the corre- 
lation coefficient between these two quantities is 0.7, while 
the correlation coefficients between the maximum flare in- 
dex and the average values of Dd and Cd are approximately 
-0.7. Interestingly, the correlation is 0.5 if Dg is averaged 
over the same 50 h time interval. 

From Table [1] we can see that there is a large scatter 
of correlation coefficients between flare indices computed 
over windows of 50 hours and ly. On the other hand, we 
noted that ly shows more rapid variations than the fractal 
and multifractal parameters. Clear and strong increases of 
its value correspond to phases of particularly frequent and 
intense activity, as we can see around 100 hours and 120 
hours from the beginning of our observation in NOAA 9236 
and NOAA 9026, respect ively (Fig. [11(a) and ( b]) 

As reported also by IVlahos fc Georgoulis I ( 2004[ ) oe 
does not show a particular trend with time, but it is in- 
teresting to notice that the average value estimated during 
the observation time for our superactive regions, except for 
NOAA 9733 and NOAA 10486, is between -1.6 and -1.8 (see 
Table [TJ last column). 



5. Discussion and conclusions 

We investigated the temporal evolution of the generalized 
fractal dimension, the contribution diversity, the dimen- 
sional diversity, the volume occupied by the potentially un- 
stable regions and the power law index of the distribution 
function of the total free energy in 14 superactive regions. 
The first three quantities provide information about com- 
plexity of the magnetic field in the photosphere, while the 
other two provide information about magnetic field organi- 
zation in the upper layers of the solar atmosphere. In par- 
ticular, we investigated the variations of these indicators 
with the emergence of magnetic flux. With the exception 
of the fractal dimension, we have also employed these indi- 
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cators for the first time for statistical studies aimed at the 
investigation of their correlation with flare productivity. 

We found that fractal dimension and multifractal pa- 
rameters are anticorrelated. During the emergence phase, 
the fractal dimension rapidly increases and multifractal pa- 
rameters decrease. The opposite was observed during the 
decay phase. In the 'mature' phase, fractal dimension stays 
fairly constant, while larger variations have been observed 
for multifractal parameters. Movies reveal that these varia- 
tions are usually due to the emergence of new magnetic flux. 
These findings are in agreement with results presented by 
iConlon et al. 1 1)200 8') , who also observed a rapid increase of 
multifractal parameters at the very beginning of the emer- 
gence phase. We do not find such an increase in our results, 
since we analyzed images taken after the observation of 
the first flare, that is, after the emergence of a large part 
of the magnetic flux. These results confirm the fact that 
magnetic flux emerges fragmented in the photosphere and 
then it evolves in more and more compact and complex re- 
gions. Moreover, the evolution of the magnetic field in the 
photosphere seems to influence the volume occupied by the 
magnetic discontinuities, that is the value of and the 
distribution functions of the energy of unstable volumes, 
only in part. Of course, the more flux that emerges, the 
greater is the probability of formation of unstable regions. 
Nevertheless, the complexity of the magnetic field in the up- 
per layers of the solar atmosphere can be l argely increased 
by ph otospheric magnetic field dynamics (|Choe fc Cheng I 
|200C1| ). that fractal and multifractal studies do not properly 
take into account for. 

We also analyzed the temporal evolution of the distri- 
bution functions of the total energy of unstable volumes. 
We found that for most of the observed regions, during the 
emergence phase the distribution function peaks as well as 
the linearity range shift toward higher energy values, while 
the slope of the linear fit (parameter as) does not vary 
significantly with the evolution of the region. 

In 6 ARs we found that the flare index is positively cor- 
related with Dg and anticorrelated with C'd and Dd- This 
indicates that in many cases the photospheric complexity is 
a strong condition for flare production. On the other hand, 
the absence of any correlation in the other ARs confirms 
that the evolution of the photospheric magnetic field inves- 
tigated by the fractal and multifractal parameters might 
not clearly reflect the real complexity of the AR in the up- 
per levels of the solar atmosphere. 

In agre ement with previous works of iMcAteer et al. I 
(|2005h and lAbrame^AoT (l2005bn . we also found that more 
com plex regions are more likely to develop flares. In partic- 
ular, lAbramenko I (|2005b[ ) found that the power spectrum 
index of ARs magnetograms and the daily flare index are 
correlated and that their variation can be described by a 
power law. We do not find such a law in our data (see FigO, 
most likely because our sample includes only superactive re- 
gions. Nevertheless, we also find good correlation between 
flare index and average fractal dimension of a region, and 
good anticorrelation between flare index and average mul- 
tifractal indicators. 

For ly we did not find a significant correlation with the 
flare index computed over time windows of 50 hours, but we 
noted a clear increase of this parameter during the phases 
when the ARs showed a burst of flares. This suggests that 
the presence of magnetic discontinuities plays an important 
role in the determination of the conditions useful for the oc- 



currence of magnetic reconnection. Moreover, the fact that 
the flares occur more frequently while the volume of po- 
tentially unstable regions increases, indicates that the flare 
itself does not provide any mechanism to stop the forma- 
tion of the instability regions. Therefore, we believe that 
the source of the current sheet formation should not be 
searched for only inside the discontinuities, but in the whole 
AR evolution. 

Although we did not observe any significant variation 
with time of as, we obtained avera ge values of this indi- 
cator below the values obtained by IVlahos fc Georgoulis I 
(2004') over two quiet ARs. This confirms the hypothesis of 
these authors according to which the steeper the slope of 
the free energy distribution function of unstable volumes of 
an AR, the higher is the probability of flare production. 

We conclude that among the complexity parameters an- 
alyzed and for the time resolution of the dataset, the gen- 
eralized fractal dimension Dd and the potential unstable 
volume of magnetic discontinuities ly are the most suitable 
measurements to investigate the flare activity over longer 
(50 hours) and shorter (few hours) time intervals, respec- 
tively. Moreover, Cd is less useful and is redundant if Dd 
is already estimated. This result is in contrast with result s 
obtained in p reviou s studies, such as Georgoulisi (|2005[ ): 



lAbramenkol (|2005af) : lAbramenko et al. I (|2008F which re- 



ported variations in the degree of multifractality of ARs 
with flare productivity. These discrepancies are likely due to 
the different spatial and temporal resolution of our dataset 
and to the e mployment of differ ent multifractal estimators. 
In particular. lGeorgoulis I (|2005f ) compared results obtained 
with different estimators and found that the structure func- 
tions and their inertial-scale scaling exponents show larger 
temporal variations than the spectrum of the generalized 
dimension, that is, the variation of Dq (defined in eqn|2]) 
with q, with flare occurrence. Since the multifractal spec- 
trum is a function of Dq, our results confirm, through an 
analysis of a larg er sample of data, results presented in 
iGeorgoulis I (j2005D . The different spatial and temporal res- 
olution of the data employed in this study can also con- 
tribute to explain the discrepancy of our results. 

A better spatial and temporal resolution is needed to 
investigate in detail our findings. Therefore, the analyses of 
MDI high resolution or HINODE/SOT data could provide 
a good test for our results. Larger samples of data and sta- 
tistical techniques more sop histicated than the cor relation, 
such as the one proposed bv lBarnes fc Lekal(|2008h . should 
be employed to evaluate the ability of complexity indicators 
to forecast flares. 
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Appendix A: Parameter evolution 

In Figs. IA.lt ?? we report the evolution of the flare activity, total 
magnetic flux ^tot, generalized fractal dimension Dg, contribution 
diversity Cj, dimensional diversity D^i and volume occupied by the 
potentially unstable regions Iv for each of the 14 superactive regions. 
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Fig. A.l. Flare activity, ^tot, Ds, Cd, Dj, and /„ vs. time for NOAA 9026 
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Fig. A.2. Flare activity, <^tot, Ds, Cd, Da and /„ vs. time for NOAA 9236 
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